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Consider a system consisting of n d-dimensional quantum 
particles (qudits), and suppose that we want to optimize the 
entanglement between each pair. One can ask the following 
basic question regarding the sharing of entanglement: what is 
the largest possible value Emax {n, d) of the minimum entan- 
glement between any two particles in the system? (Here we 
take the entanglement of formation as our measure of entan- 
glement.) For n — 3 and d = 2, that is, for a system of three 
qubits, the answer is known; Emax{3,2) — 0.550. In this 
paper we consider hrst a system of d qudits and show that 
Emax{d, d) > 1. We then consider a system of three particles, 
with three different values of d. Our results for the three- 
particle case suggest that as the dimension d increases, the 
particles can share a greater fraction of their entanglement 
capacity. 
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Quantum entanglement, as exhibited, for example, in 
the sing let state (1/V2)(|01) - |10)) of a pair of qubits, 
has been the object of much study in recent years be- 
cause of its connection with quantum communication 
and quantum computation Though entanglement 
is a kind of correlation, it is known to be fundamen- 
tally different from ordinary classical correlation. One of 
the characteristic differences is this: whereas arbitrarily 
many classical systems can be perfectly correlated with 
each other — the temperature fluctuations in ten different 
cities could in principle be exactly parallel — any entan- 
glement that may exist between two quantum particles 
seems to limit the degree to which either of the parti- 
cles can be entangled with anything else For ex- 
ample, if two qubits A and B are in the singlet state, 
then neither of them can have any entanglement with a 
third qubit C, simply because such entanglement would 
require the pair AB to be in a mixed state, whereas the 
singlet state is pure. Coffman et al. js) have generalized 
this example (still considering only qubits) by allowing 
A and B to be only partially entangled, in which case 
one finds an inequality expressing a trade-off between 
the j4i?-entanglement and the AC or SC-entanglement. 

As this sort of limitation may be a fundamental prop- 
erty of entanglement, one would like to express it more 
generally. In particular, one would like to capture quan- 
titatively the limitation on the sharing of entanglement 
among arbitrarily many particles of arbitrary dimen- 
sion. The following problem offers one approach to 
such a quantitative expression. Consider a system of n 
d-dimensional quantum particles (qudits), and suppose 
that one wants each particle to be highly entangled with 
each of the other particles. We expect that there will have 



to be compromises, since increasing the entanglement of 
any given pair will probably work against the entangle- 
ments of other pairs. It makes sense, then, to ask how 
large one can make the minimum pairwise entanglement, 
the minimum being taken over all pairs |Q . In this paper 
we address this problem, taking as our measure of entan- 
glement the entanglement of formation |^,^ , which for a 
pair of qudits ranges from zero to log2 d. For a collec- 
tion of n qudits, let us call the maximum possible value 
of the minimum pairwise entanglement Emaxin, d). This 
function, if it can be found, will give us a specific quan- 
titative bound on the degree to which entanglement can 
be shared among a number of particles. We focus in this 
paper on two special cases: n = d and n = 3. As we will 
see, our results for n — 3, combined with earlier work on 
the problem, suggest that in a well defined sense the lim- 
itation on entanglement sharing becomes less restrictive 
with increasing values of the dimension d. 

Before reviewing what is currently known about 
Emaxin, d), let us recall the definition of entanglement 
of formation. For a pure state |<I>) of a bipartite quan- 
tum system, the entanglement E(^) is defined [Q as 



EW = - Y^r^og^ 



(1) 



where the 's are the eigenvalues of the density matrix of 
either subsystem. (For a pure bipartite state the density 
matrices of both subsystems necessarily have the same 
eigenvalues.) A mixed state p can always be written in 
many different ways as a probabilistic mixture of distinct 
but not necessarily orthogonal pure states: 



P 



(2) 



The entanglement of formation of p is defined ||,^ as the 
average entanglement of the pure states of the decompo- 
sition, minimized over all possible decompositions: 



Ef{p)^Mj2pjmi)- 
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As we have mentioned above, the entanglement of forma- 
tion between a pair of qudits ranges from zero to log2 d. 
Let us refer to the maximum value log2 d as the entan- 
glement capacity of a pair of qudits. 

For a pair of qubits, there is an explicit formula for the 
entanglement of formation of an arbitrary mixed state [|| . 
It is given in terms of another measure of entanglement 
called the concurrence MbI, which at this point has a 
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standard definition only for qubits.|^ In terms of the con- 
currence C, which ranges from zero to one, the entangle- 
ment of formation of a pair of qubits is Ef{p) — £{C{p)), 
where the function £ is defined by 



S{C) = h 



(4) 



h being the binary entropy function h{x) = —x log2 x — 
(1 — a;)log2(l — x). Note that £{C) is a monotonically 
increasing function, with £{Q) = and £{1) = 1. Wc will 
not be focusing particularly on qubits in this paper, but 
Eq. (^) will be useful both for summarizing previous work 
on the problem and, in a different context, for presenting 
our own results. 

We now list the results that have been obtained so far 
regarding Emax{n, d). 

1- d) — log2 d. This equation simply says 

that if there are only two particles, they can satu- 
rate their entanglement capacity; they do not have 
to share the entanglement with other particles. 

2. Ernaxi3,2) = £{2/3) = 0.550. Diir et al. g ob- 
tained this result by proving that the optimal pair- 
wise entanglement for a system of three qubits is 
achieved in the state (l/^/I)(|100> + |010> + |001)). 

3. Emaxin,2) > £{2/n). Koashi et al. [|o| showed 
that for a system of n qubits, if the state is such 
that the density matrix of each pair of particles 
is the same, then the maximum pairwise concur- 
rence is 2/n. It is conceivable (though it seems un- 
likely) that by removing the symmetry constraint 
one might be able to achieve a greater pairwise en- 
tanglement; therefore we write this result as an in- 
equality rather than an equality. 

In this paper we add two new items to the above list: 
(i) For n = d, that is, for a system of d qudits, we find 
for each value of d a specific state in which each pair of 
particles has exactly 1 "ebit" of entanglement between 
them. (For d — 2, our state reduces to the singlet state 
of a pair of qubits.) This will show that Emax{d, d) is at 
least 1 for all values of d. (ii) For n — 3, that is, for a 
system of three particles, we add to the known result for 
qubits (d = 2) and to our own result for qutrits (d = 3 



^The concurrence of a pure state of two qubits is simply 
2y/det pA, where pA is the reduced density matrix of one of 
the qubits. The concurrence of a mixed state p of two qubits is 
given by max{0, Ai — A2 — A3 — A4}, where Ai > A2 > A3 > A4 
are the square roots of the eigenvalues of p{Gy®ay)p* {(Ty^ay), 
p* being the complex conjugate of p in the standard basis and 
Oy being the usual Pauli matrix [^,^. We recall these formulas 
here for the sake of completeness but will not need them in 
the present paper. 



in item (i)) a third example with d = 7. Our results for 
the three-particle case suggest that as d increases, the 
particles can share not just more entanglement, but a 
greater fraction of their entanglement capacity. 

A system of d qudits 

Before writing down our special state of d qudits with 
arbitrary d, we illustrate our construction in the special 
case of three qutrits. Let the particles be called A, B 
and C, and let the indices i, j and k label the elements 
of orthogonal bases for the three particles, each index 
taking the values 0, 1 and 2. Our special state for this 
system is 
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where eijk is antisymmetric under interchange of any two 
indices and eoi2 = 1- This is the singlet state of three 
qutrits with respect to the group SU{S); i.e., it is the 
unique three-qutrit state (up to an overall phase factor) 
that is invariant under arbitrary transformations of the 
form U ®U ®U where U G SU{3). The density matrix 
|^)(^| is symmetric under interchange of any two parti- 
cles, so that each pair of particles is equally entangled. 
To find the pairwise entanglement, we write down the 
reduced density matrix of any pair; for definiteness we 
choose the first two particles, A and B: 



g ^ijk^i'j'k — ^{^ii'Sjji SijiSji'), 



(6) 



d being the Kronecker delta. Alternatively, we can write 
p^^ without indices as 



/« = i(/-n 



(7) 



where / is the identity operator and F is the operator 
that interchanges particles A and B: = Ni)0*l- 
The two-qutrit state p^^ is an example of a Werner state, 
that is, a state that is invariant under all transformations 
of the form U ®U where U is unitary. Werner states can 
be defined for any dxd system, and one can show [lT| that 
in any dimension the Werner states are precisely those 
states that can be written as p = al -\- bF, a and b being 
real numbers and F being defined as above. VoUbrecht 
and Werner have shown that the entanglement of for- 
mation of any Werner state is given by Ef{p) — £{c{p)), 
where c{p) = — Tr pF and £ is the function defined in 
Eq. (jj). (When c{p) is non-negative, it plays the role of a 
concurrence for Werner states.) In our case, c{p^^) = 1, 
so that the entanglement is Ef{p^^) = £{1) = 1. Thus 
each pair of qutrits has exactly one ebit of entanglement. 
This value is, by the way, the maximum possible entan- 
glement of any Werner state. 

It is a simple matter to generalize the above construc- 
tion to a system of d qudits. In that case, we have d 
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indices ii, 12, ■ • • , id, each taking values from to d — 1. 
Our special state for this system is the SU{d) singlet 
stateQ 



where 



is completely antisymmetric and eq. 



.d-l 



1. 



One can show directly that the reduced density matrix 
of each pair of particles is again a Werner state: 



P 



AB 



1 



d{d^l) 



(I-F), 



(9) 



and that the entanglement of formation of this state is 
one ebit. We thus conclude that E„iax{d,d) > 1. We 
write an inequality here simply because our state |^) may 
not optimize the pairwise entanglement; one might be 
able to do better. However, having put some effort into 
looking for better states with d = 3, we regard it as likely 
that our state is optimal in that case. 

A three-particle system 

It is interesting to compare our result for three qutrits 
with the previously studied example of three qubits [Q. 
For a triple of qubits Emax is 0.550, and we have just 
seen that for a triple of qutrits, Emax is at least 1. How- 
ever, a straightforward comparison of these numbers is 
not particularly illuminating, because qubits and qutrits 
have different entanglement capacities. We can perhaps 
make a fairer comparison by considering the ratio of 
Emax to the relevant entanglement capacity. For qubits, 
this ratio is 0.550/ log2 2 = 0.550, whereas for qutrits it 
is l/log2 3 = 0.631. Thus by this measure, qutrits are 
better able to share entanglement than qubits: they can 
share a greater fraction of their entanglement capacity. 
It is interesting to ask whether this trend will continue 
for larger values of d. That is, will (3, d)/ log2 d 

continue to increase with increasing d? 

To address this question, we consider one further case 
with three particles, namely, the case d = 7. We choose 
the value 7 because it allows us to construct a reasonably 
simple and symmetric state that exhibits large pairwise 
entanglement. In fact we consider a one-parameter family 
of states, having the following form: 
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Ak) 



(10) 



This state has been used recently by Hillery and Buzek in 
a scheme designed to probe a quantum gate that realizes an 
unknown unitary transformation lia]. 



Here Q is the set {1, 2, 4}, and all the arithmetic shown in 
the ket labels is mod 7, the basis states of each particle 
being labeled by the integers 0, . . . , 6. We take a and 
h to be real and positive, with + 36^ = 1 to ensure 
normalization. Thus the state \Q is completely specified 
once the value of a is given. 

We have chosen Q to consist of the quadratic residues 
mod 7, that is, the elements of {1, 2, 3, 4, 5, 6} that can be 
written as mod 7 for some integer x. The properties 
of quadratic residues JTst tend to minimize the overlap, 
in each particle's state space, between terms in Eq. ( p^ ) 
with different values of j. (For this it is important that 
7 is a prime of the form AN — 1 with integral N .) 

A particular symmetry of \C,) shows immediately that 
the pairs AB^ BC, and CA are all equally entangled. If 
we define a new summation index A;' by A; = 2k' mod 7, 
then Eq. (nw becomes 



1 

(11) 



fc'ec 



where we have used the invariance of Q under multipli- 
cation by 2 mod 7. But Eq. (|ll]) differs from Eq. (^ in 
that the ket labels have been cyclically permuted. Thus 
\Q is invariant under a cyclic permutation of the parti- 
cles, and it follows that each pair is equally entangled. 

To write down the density matrix of one of the pairs, 
say BC, it is helpful to re-express Eq. ( p^ ) in yet another 
form, changing the index j in the jfc-sum to j' = j + k 
and then relabeling j' as j: 

1 ^ 

10 = 7f E + ^ E 1^' ^ + + 3fc>) • (12) 



J=0 



keQ 



The density matrix of BC is the trace of |C)(CI over par- 
ticle A, which we can write as 



1 ^ 

PBC = y E I^J')^' 
j=0 



where |sj), defined by 



\sj) = a\j,j) +bY,\j + k,j + 3fc), 
keQ 



(13) 



(14) 



is the state of BC associated with the state \ j) of A. 

In order to find the entanglement of formation E{pBc), 
we need to consider pure-state decompositions of psc 
and find their average entanglements. Now, any pure 
state 1/3) in such a decomposition must be a linear com- 
bination of the seven orthogonal states \sj) that make 
up pbc\ that is, it must lie in the seven-dimensional sub- 
space Ti spanned by {|sj)}: 



-E/3: 



(15) 
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where J^j IPjl"^ = 1- The problem of finding E{pBc) is 
simphfied by two facts: (i) E{pBc) cannot be smaller 
than the smallest entanglement of any S H; that is, 
E{pBc) > min^ E{(i). (ii) Given any state \(3) G TL, one 
can generate an entire decomposition of psc in which 
every element has the same entanglement as \(3). (We 
prove this assertion in the following paragraph.) To- 
gether, these two facts imply that the entanglement of 
formation of p^c is equal to m\np E{(3). Thus it is suf- 
ficient to find a single minimally entangled pure state in 
the subspace 7i occupied by psc- 

To generate a decomposition of psc from a given state 
1/3) G Ti, we apply a set of local unitary transformations 
to 1/3) ; such transformations are guaranteed not to change 
the entanglement. We start by defining two basic single- 
particle transformations S and T: 

s\]) = o^\])\ m^\] + i), (16) 

where w — e^'^*/'' and, as always, the addition in the ket 
label is mod 7. In terms of these basic operations, we 
define a pair of two-particle transformations U and V: 



One can show that U\sj) 
from which it follows that 



V = T 



E>T. 

and V\sj) 



(17) 



41 ELo E'=o VPU-^m (/3| C/-"l/-P 

= l^jPEp \Sj+p){Sj+p\ = PBC- 



(18) 
(19) 
(20) 



We have thus produced the desired decomposition of 

PBC- 

It remains, then, to find the smallest possible value of 
E{j3). For the special case where a = b ^ 1/2, each state 
\sj) has exactly two bits of entanglement, but it happens 
that certain linear combinations of the states \sj) have 
slightly smaller entanglement. Using numerical mini- 
mization, we find that for this case, mini?(/3) — 1.9933. 

Of course we are free to choose the value of 
a as we please, and it turns out that we maxi- 
mize the entanglement of formation by choosing a = 
0.461, in which case b — 0.512. For this value 
of a, we find numerically that the minimum E{(3) 
is 1.9944, obtained both for the simple case |/3) — 
\sj) and for certain non-trivial linear combinations. 
[One such combination has coefficients Pj equal to 
(0.120, 0.197, 0.689, 0.259, -0.468, -0.275, -0.332), and 
the others we have found are all related to this one by 
permutations and phase changes.] We conclude, then, 
that for the state |C) with a = 0.461, the entanglement 
of formation between each pair of particles is 1.9944, and 
therefore E^ax{5,7) > 1.9944. 

This result gives us another data point as we con- 
sider the dependence of the ratio iJmaa; (3, d)/ logj c? on 
the dimension d. The following table summarizes what 
we know so far about the case n — 3. (For d — 3 and 
d = 7, the values given are lower bounds.) 



d 


^max l,*-*, ^} 


i?L"^(3,d)/log2d 


2 


0.550 


0.550 


3 


1.000 


0.631 


7 


1.994 


0.710 



In the limit as d goes to infinity, we wonder what value, 
if any, the ratio Ernax{i, c^)/log2 d approaches. It is con- 
ceivable that the limit is 1, but it is equally conceivable 
that it is some smaller constant. Either answer would be 
interesting. If the limit of Emax{n, o?)/ log2 d is 1 for all 
values of n, then one could reasonably say that entangle- 
ment can be shared freely in an infinite dimensional state 
space. 

We note that although in this paper we have focused 
on the entanglement of formation, there exist other sound 
measures of entanglement, and it is surely a good idea, 
in trying to quantify the restrictions on the sharing of 
entanglement, to keep in mind alternatives such as the 
relative entropy of entanglement [^ and the generalized 
concurrence of Rungta et al. [^ . At the present stage of 
investigation, it is not clear which measure or measures 
will yield the most elegant quantitative expressions of the 
limitations on entanglement sharing. 
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